Broadband dielectric response of glycerol 
and propylene carbonate: a comparison 

P. Lunkenheimer, U. Schneider, R. Brand, and A. Loidl 

Experimentalphysik V, Universitat Augsburg, D-86135 Augsburg, Germany 



Abstract. Dielectric data on glycerol and propylene carbonate covering 18 decades of 
frequency are presented and compared to each other. Both materials exhibit qualita- 
tively similar behavior except for marked differences in the high-frequency region just 
below the boson peak. The results on both materials are consistent with the mode 
coupling theory of the glass transition. 



I INTRODUCTION 

In recent years, a variety of new theoretical and phenomenological approaches 
of the glass transition (e.g., [1-3]) stimulated new experimental investigations es- 
pecially of the high-frequency dynamics of glass-forming liquids [4-15]. Maybe the 
most controversially discussed theoretical approach of the glass transition is the 
mode coupling theory (MCT) [1], which explains the glass transition in terms of a 
dynamic phase transition at a critical temperature T c significantly above the glass 
temperature T g . For frequencies in the GHz-THz region, MCT predicts an addi- 
tional contribution, now commonly termed fast /3-relaxation. This frequency region 
was mainly investigated by neutron and light scattering experiments [5-8] . But re- 
cently, by combining various techniques, our group was able to obtain continuous 
dielectric spectra on glass-forming liquids extending well into the relevant region 
[10-15]. For glass-forming glycerol and propylene carbonate, spectra covering 18 
decades of frequency and extending well into the THz range were obtained [14,15]. 
Glycerol is a rather strong [16] hydrogen-bonded glass-former with a fragility pa- 
rameter of m ~ 53 [17]. In contrast, propylene carbonate (PC) can be characterized 
as a fragile (m ~ 104, [17]) van der Waals liquid. In the present paper we present 
a comparison of the results on both materials. For the experimental details the 
reader is referred to our earlier publications [11,14,15]. 



II RESULTS AND DISCUSSION 



Figure 1 shows the e"(z/)-spectra for glycerol and PC in the whole accessible 
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FIGURE 1. e"(v) of PC (a) and glycerol (b) at various temperatures. The dashed lines are fits 
of the a-peak with the CD function. The solid lines are fits with a phenomenological ansatz (see 
text). The dash-dotted lines indicate a power law increase at v < 1 THz. The FIR results in (a) 
have been connected by a dotted line to guide the eye. 

frequency range. e"(u) exhibits the typical asymmetrically shaped a-relaxation 
peaks shifting through the frequency window with temperature. The data agree 
with the results of earlier dielectric investigations [4,9,18-20] which were restricted 
to smaller frequency and temperature ranges. The dashed lines in Fig. 1 are fits of 
the a- relaxation region with the empirical Cole-Davidson function [20] . A good fit 
of the peak region was achieved. The temperature dependence of the parameters of 
the CD function, < tqd >= tcdPcd, Pcd, and Ae = e s —eoo [20] is shown in Fig. 2. 
In most respects, the results are in accord with previously published data [4,18-21] 
extending them to lower or higher temperatures. v T (T) can be parameterized using 
the Vogel-Fulcher-Tamman (VFT) equation, < t cd >= T exp[DT VF /(T - T VF )\ 
with Typ = 132 K, D = 6.6 for PC and T VF = 131 K, D = 16.5 for glycerol [dashed 
lines in Figs. 2(a) and (b)]. The values of the strength parameter D characterize PC 
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FIGURE 2. Parameters from simultaneous fits of e'(u) and e"(y) with the CD function. The 
dashed and solid lines in (a) and (b) are fits with the VFT law and the MCT critical behavior, 
respectively. The lines in (b) and (c) demonstrate a possible saturation of Pcd(T) at high 
temperatures. The lines in (e) and (f) are calculated according to MCT (see text). 

as fragile and glycerol as intermediate in Angell's classification scheme [16]. At high 
temperatures small deviations from VFT behavior show up, similar to those seen in 
earlier work [4,19,21]. The simplest version of MCT, the idealized MCT [1], predicts 
a critical behavior of the a-relaxation timescale, v T = l/(2n < tqd >) ~ (T — T c ) 7 
with 7 = l/(2a) + 1/(26). Here a and b are the power law exponents of the s"{v)- 
minimum (see below). The solid lines in Figs. 2(a) and (b) are fits with the MCT 
prediction with 7 = 2.72, T c = 187 K for PC and 7 = 2.33, T c = 262 K for glycerol, 
fixed to the values obtained from the analysis of the e"(z/)-minimum (see below). 
Only the data above 200 K for PC and above 280 K for glycerol have been used 
for the fits which led to a good agreement in this temperature region. For lower 
temperatures deviations show up. This can be understood within the extended 
MCT [1] where the structural arrest found in idealized MCT for T < T c is overcome 
by the introduction of thermally activated hopping processes. 

For both materials, the width parameter (3cd{T) [Figs. 2(c) and (d)] increases 
linearly below about 200 K for PC and 300 K for glycerol. For T > T c , MCT 
predicts a temperature independent spectral form of the a-process. Indeed, a ten- 
dency to saturate is seen for (3 C d{T) in Figs. 2(c) and (d), however at temperatures 
somewhat above T c , only. But a reasonable description of the a-peaks for T > T c is 
also possible with a constant 0cd = 0.63 for glycerol [13] and (3 C d = 0.80 for PC. 




The Ae(T) of PC and glycerol [Figs. 2(e) and (f)] both decrease with temperature. 
According to MCT, Ae should be temperature independent for T > T c . In addition, 
for T < T c , Ae = ci+ c 2 [(T c - T)/T^ follows from extended MCT. In Figs. 2(e) 
and (f) the lines were calculated using the MCT prediction. The agreement of fits 
and data is rather poor, but due to the large experimental uncertainties no definite 
conclusion can be drawn. At least it seems that Ae assumes a weaker temperature 
dependence above T c . 

At frequencies about 2-3 decades above the peak frequency, for both materials 
an excess wing shows up as power law, e" ~ v~ h with b < [3cd (Fig. 1). The 
exponent b decreases with decreasing temperature as found previously for various 
glass-forming materials [9]. For PC, at 203 K and for glycerol at 272 K, the excess 
wing has merged with the a-peak. The excess wing seems to be a universal feature 
of glass-forming liquids [3] but, up to now, no consensus concerning its microscopic 
origin was achieved. For a more detailed treatment of the excess wing in glycerol 
and PC, see [15]. 

At frequencies in the GHz-THz range, e"(u) exhibits a minimum for both ma- 
terials [10,11]. A magnified view is given in Fig. 3. Clearly this minimum is too 
shallow to be described by a simple superposition of a-peak (or excess wing) and 
boson peak [10,11]. Near 1 THz, for both materials the FIR measurements reveal 
a peak or shoulder, that can be ascribed to the boson peak, well known from light 
and neutron scattering experiments [7,8]. Both materials exhibit a significantly 
different behavior in the region between minimum and boson peak: For PC only 
one power law, e" ~ u a , is seen at frequencies v > z/ min forming simultaneously the 
high frequency wing of the ^''-minimum and the low frequency wing of the boson 
peak. For T — > T g , a linear behavior is approached as indicated by the dash-dotted 
lines in Figs. 1(a) and 3(a). In marked contrast, for glycerol two regimes can be 
distinguished: Just above f m in, there is a rather shallow increase of e"{v) but at 
higher frequencies a very steep increase appears, approaching e" ~ u 3 for T — > T g 
[dash-dotted lines in Figs. 1(b) and 3(b)]. 

For both materials, the e"(z/)-spectra between a- and boson peak can be described 
by the phenomenological ansatz cpv^ 13 + q>z/~ & + e c + c^v ' 3 + c n v n (solid lines in Fig. 
1). The first term takes account of the high-frequency wing of the a-peak with the 
parameters chosen to achieve a smooth transition to the CD-fits. The second term 
describes the excess wing. In order to explain the shallow minimum, a constant 
loss term, e c , as proposed by Wong and Angel [22] and a e" ~ z/ a3 increase were 
assumed. Both contributions seem to be universal features in the high-frequency 
response of glassy ionic conductors [23]. Finally, the v n power law, with n > 1, 
takes account of the low frequency wing of the boson peak. In this way, quite 
impressive fits over 17 decades of frequency are possible. Of course this could be 
expected considering the large amount of parameters, but nevertheless it cannot 
be fully excluded that such a simple superposition ansatz is correct. 

MCT predicts for the minimum region significant contributions from the fast 
/9-process. Within idealized MCT, above T c , the minimum region can be approxi- 
mated by the interpolation formula [1] e"{v) = e m i a [a(v> I V min )~ b +5(z/ /u mia ) a ]/ (a+b). 
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FIGURE 3. e"(v) of PC (a) and glycerol (b) at high frequencies. The solid lines are fits with 
the MCT prediction (see text). The dashed lines demonstrate v a behavior for T < T c . The 
dash-dotted lines indicate a power law increase at v < 1 THz. The FIR results in (a) have been 
connected by a dotted line to guide the eye. 



Here v min and e" min denote frequency and amplitude of the minimum, respec- 
tively. The exponents a and b are temperature independent and constrained by 
A = T 2 (l - a) /r(l - 2a) = T 2 (l + b) /T(l + 2b) where T denotes the Gamma func- 
tion. This implies a < 0.4, i.e. a significantly sublinear increase of e"(u) at v > v min 
is predicted. For PC, a consistent description of the ^''-minima at T > 193 K is 
possible with A = 0.76 (a = 0.29, b = 0.5) [11,15] [solid lines in Fig. 3(a)]. The ob- 
tained A is consistent with the results from other measurement techniques [8,24,25]. 
The fits provide a good description of the data up to the boson-peak frequency. 
Quite a different behavior is seen in glycerol [solid lines in Fig. 3(b)]: Here the 
MCT fits (A = 0.705, a = 0.325, b = 0.63) [10,11] are limited at high frequencies 
by the additional steeper increase, mentioned above. It was argued [10,11], that 
this deviation in glycerol is due to the vibrational boson peak contribution which 
is not taken into account by the MCT interpolation formula. In PC this contribu- 
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FIGURE 4. Temperature dependence of £^ in , t-Wn, and v T plotted in representations that result 
in straight lines according to the predictions of MCT. The solid lines extrapolate to a T c of 187 
K for PC and 262 K for glycerol. 



tion seems to be of less importance. This is in accord with the finding of Sokolov 
et al. [26] that the amplitude ratio of boson peak and fast process is largest for 
strong glass formers, glycerol being much stronger than PC. However, it has to be 
noted, that in the present results, at least near T g this ratio is quite similar for 
both materials. 

For T > T c , idealized MCT predicts the following critical temperature depen- 
dences: v mm ~ (T - T c ) l K 2a \ e" mm ~ (T - T c )h, and v T ~ (T - T c )\ In Figure 4 
representations have been chosen that should lead to straight lines that extrapolate 
to T c . Indeed, all three data sets can be described consistently with a T c ps 187 
K for PC [11,15] and T c 262 K for glycerol [10] as indicated by the solid lines. 
These values lie in the same range as the T c -value obtained from other techniques 
[7,8,24,25,27]. For temperatures near T c the data partly deviate from the predicted 
behavior. Within MCT this can be ascribed to a smearing out of the critical be- 
havior near T c due to hopping processes which are considered in extended versions 
of MCT only [1]. In addition, the above critical temperature dependences of s" min , 



Vmini and v T should be valid only for temperatures not too far above T c . Therefore 
the proper choice of the temperature range to be used for the determination of T c 
is difficult which may lead to some uncertainties concerning the value of T c . 

For T < T c , MCT predicts the occurrence of a so-called "knee" at v — showing 
up as a change of power law from e" ~ u a at v > to e" ~ v at v < Uk [1]. In 
Fig. 3(a) the dashed lines suggest a e" ~ v a behavior for frequencies just below 
the boson peak. For lower frequencies the experimental data exhibit a downward 
deviation from the z/ a -lines. This may be indicative of the "knee" predicted by 
MCT, but, clearly, the data are not of sufficient quality to demonstrate conclusively 
the existence of the "knee". 

Ill CONCLUSIONS 

Dielectric data on glass-forming PC and glycerol in an exceptionally broad fre- 
quency range have been presented and compared to each other and with the pre- 
dictions of MCT. Concerning the a-relaxation and the excess wing, both materials 
behave qualitatively similar. The parameters describing the spectral form and dy- 
namics of the a-relaxation are in rough agreement with the MCT predictions but 
clearly, Fig. 2 is not very convincing in this respect. The high-frequency response 
of both materials provides clear evidence for the presence of fast processes in the 
GHz-THz region. For glycerol an additional very steep increase towards the boson 
peak seems to be superimposed to the shallow e"{y) -minimum. For the more fragile 
PC a smooth sublinear increase of s"{y) is seen up to the boson peak frequency. For 
both materials the frequency and temperature dependence in the minimum region 
can be consistently described by MCT. Admittedly, also alternative explanations 
may be possible. But even in its idealized form, MCT provides a consistent picture 
for a large variety of experimental facts. Obviously, more sophisticated approaches 
within MCT, e.g. including hopping processes [1], the boson peak [28], or orienta- 
tional degrees of freedom [29] are necessary to eliminate deviations still seen in the 
present work. 
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